On the non-existence for quantum LDPC codes of 
type IEEE802.16e with rates 1/2 and 2/3B 
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I. Introduction 

In this paper, we discuss a construction of CSS codes 
derived from pairs of practical irregular LDPC codes. Inter- 
section studies between quantum error-correcting codes and 
LDPC codes are current hot research topics [1], [3], [4], [5], 
[6], [7], [8], [9]. One of aims of quantum error-correcting 
code theory is to construct quantum error-correcting codes 
with small length and high error-correcting performance. The 
reason why LDPC codes are chosen is their high error- 
correcting performances as classical error-correcting codes. 
In fact, LDPC codes have almost achieved a theoretical 
error-correction limit, called a Shannon limit, with various 
conditions, for example, code rate, code length, communicatin 
channel, and others. This fact induces naturally derived the 
quantum error-correcting codes. 

In the researches, we encounter a problems to hold a 
condition, which we call it the "twisted condition", for a pair 
of classical codes to derive a quantum code. The condition 
requires a geometrically orderly structure to two classical 
code spaces. More explicitly, the dual code of one of two 
codes must be contained in the other code. On the other 
hand, a complicated structure, i.e. a random structure, is 
suitable for achieving high error-correcting performance. The 
construction theory of classical codes for a CSS code is to 
satisfy conflicting, orderly and complicated, requirements. We 
re-write the two conditions in terms of parity-check matrices. 
The twisted condition implies the orthogonality of their parity- 
check matrices. A complex structure implies that non-zero 
elements in the parity-check matrices should be irregularly 
arranged. 

The previous researches on quantum LDPC codes have 
focused on quantum regular LDPC codes [1], [3], [7], [8] 
as same as that of classical LDPC codes. But now, a main 
stream of classical LDPC code researches is a construction 
method for irregular LDPC codes. The crucial reason why 
the main research stream has changed is error-correcting 
performance. In fact, the performances of irregular LDPC 
codes are extremely better than those of regular LDPC codes. 
This fact is theoretically proved by density-evolution theory 
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and the Gaussian approximation method [10], [11]. 

Our design of irregular LDPC codes is based the design 
written in the standardization of IEEE802.16e. It is need- 
less to say that irregular LDPC codes which are chosen 
in IEEE802.16e are designed well-considerably and show 
high error-correcting performance with practical length. Our 
research has tried to make a CSS code with a pair of LDPC 
codes of type IEEE802.16e. To our regret, we proved that it 
was impossible to construct a CSS code if one of classical 
codes was of type IEEE802.16e with rate 1/2 and 2/3B. We 
would like to report the discussion on its impossibility in this 
paper This is the first paper to analyze the possibility of a 
CSS code construction by using two irregular LDPC codes 
which are practically useful. 



11. Preliminaries 



A. CSS codes 



Denote the binary field by F2. CSS (Calderbank-Shor- 
Steane) codes are quantum codes constructed by a pair of 
classical linear codes C and D which satisfy the following 
condition (T), called twisted relation: 

(T) C C, 

where is the dual code of D. The dual code of D is 
defined by the following: 

:= {d'\d X d'^ = 0,Vd G D}, 

where d'^ is the transposed vector of d'. A function (, ), 
defined by {x,y) :— x x , gives an inner product over the 
binary field F2. The dual code is regarded as an orthogonal 
space of D via (, ). Note that is equal to a linear code 
generated by a parity-check matrix Ho of D. In other words, 

= {xHd\x e F^"}, 

where m is the column size of Hd. 

Details of the construction of a "quantum" CSS code by 
such a pair of "classical" codes C and D is written in [13]. 

Remark 2.1 : In general, a pair of linear codes forming a 
CSS code is denoted by Ci and C^. For simplicity, we use 
C (resp. D) instead of Ci (resp. C^). 



B. Quasi-Cyclic LDPC Codes 

A quasi-cyclic LDPC code with circulant matrices (QC- 
LDPC code) is defined by a binary matrix H of size m-by-n, 
where n is the length of the code and m is the number of 
parity check bits in the code [2]. The number of systematic 
bits is k = n — m. The matrix H is defined as: 



H = 



J-1,0 



PlA 



Po,L-l 
Pl,L-l 

Pj-1,L-1 



where Pij is one of a set of z-hy-z circulant permutation 
matrices or a z-by-z zero matrix and z is a positive integer. 
The matrix H is expanded from a binary base matrix Hh of 
size J-by-L, where m = zJ and n = zL. The base matrix 
is expanded by replacing each 1 in the base matrix with a 
^-by-2; circulant permutation matrix, and each with a ^-by- 
z zero matrix. The circulant permutations used are circular 
right shifts, and the set of permutation matrices contains 
the .2X2; identity matrix and circular right shifted versions 
of the identity matrix. Because each permutation matrix is 
specified by a single circular right shift, the binary base matrix 
information and permutation replacement information can be 
combined into a single compact model matrix Ti. The model 
matrix H is the same size as the binary base matrix Hb, 
with each binary entry of the base matrix Hb replaced 
to create the model matrix Ti. Each in Hb is replaced by 
a symbol oo to denote a z x z all-zero matrix, and each 1 
in Hb is replaced by a circular shift size p{i,j) > 0. The 
model matrix H can then be directly expanded to H. Denote 
a circulant permutation matrix of single shift by /(I), i.e. 



7(1) = 



V 1 



1 



and I{b) = /(l)**. A circulant permutation Pj_j in is 

mhj)). 

Example 2.2: Let be a parity-check matrix of size 6-by- 
10 with the following form: 



H 



Then H defines a quasi-cyclic LDPC code with circulant 
matrices of size 2. The base matrix Hb associated to H is 
a matrix of size 3-by-5: 
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Denote a set {0, 1, . 
matrix is a matrix over \z] 



The model matrix H associated to i7 is a matrix of size 3-by- 
5: 

1 00 

1 oo 
,z — 1} U {oo} by [z]ac- A model 
3. We define a operator — on [z]oo- 
For a,b G {0, 1, . . . , 2 — 1}, we define a — 6 as an usual 
integer operation modulo z. For a G {0, 1, . . . , 2; — 1}, we 
define a — 00 — 00 — a = oo — 00 = 00. For vectors v 

. vi, . . . , vl-i) and u = {uq, ui, . . . , ul-i), we define v — 
u = {vq — uo,vi — ui,. . . ,vl-i — ul-i). For example, we 
have (00, 0, 1, 0, 00) — (1, 0, 1, 00, 0) = (00, 0, 0, 00, 00). 

C. Quasi-Cyclic LDPC codes of type IEEE802.16e 

In the standardization IEEE802.16e, quasi-cyclic LDPC 
codes are chosen as one of optimal error-correcting codes and 
six model matrices are written to define quasi-cycUc LDPC 
codes. For each model matrix, 19 kinds of the size of cireculent 
matrices are chosen. Totaly, 1 14 quasi-cycUc LDPC codes are 
chosen in IEEE802.16e standards. These six model matrices 
are designed to satisfy the following conditions. In this paper, 
we call a quasi-cyclic LDPC code (a LDPC matrix, a base 
matrix, and a model matrix) which satisfies the following 
condition of type IEEE802.16e. 

Remark 2.3: The conditions below is not a characterization 
but a generalization of the 6 model matrices in IEEE802.16e. 
In fact, the set of model matrices which satisfy the condition 
below contains the 6 model matrices in IEEE. 

The base matrix Hb of a quasi-cyclic LDPC matrix is 
partitioned into two sections, where Hbi corresponds to the 
systematic bits and Hb2 corresponds to the parity-check bits, 
such that Hb = [(-fffci)jx(L-j) l(^^&2)jx j]- 

Section Hb2 is further partitioned into two sections, where 
hb is a vector whose weight is three, and iJ^2 ^ dual- 
diagonal structure with matrix elements at row i, column j 
equal to 1 for i = j, 1 for i = j -\- 1, and elsewhere: 



Hb2 = [hblH',^] 



hb{0) 
hb{l) 



_ hbiJ-1) 



1 

1 1 

1 



1 1 
1 



The base matrix has hb{0) = l,hb{J — 1) = 1, and a third 
value hb{j) = 1 with some < j < J — 1. In particular, 
the non-zero sub-matrices are circularly right shifted by a 
particular circular shift value. Each 1 in ff^2 is assigned a 
shift size of 0, and is replaced by a 2 x 2; identity matrix when 
expanding to H. The two located at the top and the bottom of 
hb are assigned equal shift sizes, and the third 1 in the middle 
of hb is given an unpaired shift size. 

Example 2.4: We pick up two base matrices Hb(i/2) and 
Hb{2/3B) written in the standardization IEEE802.16e. The 



matrix H;,(i/2) defines quasi-cyclic LDPC codes of rate 1/2 
and has the following form: 

" 011000001100110000000000 " 
010001110001011000000000 
000111010001001100000000 
101000001100000110000000 
001000100110000011000000 
_ 000011010001100001100000 

^6(1/2) - 001100000110000000110000 ■ 

011000100100000000011000 
100011010001000000001100 
000001010011000000000110 
001 100001 10000000000001 1 
100001010001100000000001 

The matrix -ffb(2/3S) defines quasi-cyclic LDPC codes of rate 
2/3 and has the following form: 

" 101010101010101011000000 " 
010101010101010101100000 
101010101010101000110000 
_ 010101010101010100011000 

J^b{2/3B) - 101010101010101000001100 ■ 

010101010101010100000110 
101010101010101010000011 
010101010101010110000001 

D. Preliminaries of the twisted condition for Quasi-Cyclic 
LDPC codes 

Recall that a CSS code is defined by two classical linear 
codes. Denote the two classical codes by C and D. We denote 
a quasi-cyclic LDPC matrix, the model matrix, and the base 
matrix associated with C by He, ^^6(c)' ™d Tic, respectively. 
Similarly, we use the notation Hn, and Tio to denote 

a quasi-cyclic LDPC matrix, the model matrix, and the base 
matrix associated with D, respectively. 

In [7], a necessary and sufficient condition for QC-LDPC 
codes C and D to satisfy (T) in terms of the model matrices 
has been obtained. We quote the necessary and sufficient 
condition from [7]. Denote j th rows of the model matrices 
C and V by Cj and dj , respectively. We call a vector v = 
{vq, vi, . . . , over {0, 1, . . . , z — 1} U {oo} multiplicity- 

even if each symbol except for oo appears even times in 
{vo, wi, . . . , vl-i}- For example, a vector (0, 0, 1, 1, 1, oo, 1) 
is multiplicity-even. 

Theorem 2.5 ([7] Prop. 3.1. and Sec. IV(D)): Let C and 
D be quasi-cyclic LDPC codes. The codes C and D satisfy 
(T) if and only if Cj — dk is multiplicity even for any row Cj 
of the model matrix He and any row dk of the model matrix 
Hd- 

III. A Necessary Condition for the Twisted 
Condition in Terms of the Base Matrices 

Proposition 3.1: Let C and D be quasi-cyclic LDPC codes 
with circulant matrices of size z. Let y be a positive integer 
such that y divides z. Let C (resp. D') be a quasi-cyclic 
LDPC code with the model matrix same as C (resp. D) and 



with circulant matrices of size y. If C and D satisfy the twisted 
condition then C" and D' satisfy the twisted condition. 

Proof: For any positive integers a,b, if a = b (mod z) 
then a mod b (mod y). It implies that if a vector v over [z]oo 
is multiplicity-even then the vector v is multipUcity-even as a 
vector over [y]oc- 

Let Ci (resp. di) be the ith row of the model matrix He 
(resp. Hd) of C (resp. D). By Theorem 12.51 cj — dk is 
multiplicity-even as a vector over [z]oo- Thus Cj — dk is 
multiplicity-even as a vector over [y]oo- ■ 

Corollary 3.2: Let C and D be quasi-cyclic LDPC codes 
with circulant matrices. Let H^c) (resp. Hii(^d)) be the base 
matrix of C (resp. D). Let C" (resp. D') be a linear code 
with a parity-check matrix Hb(c) (resp. H^ji^d))- If C and D 
satisfy the twisted condition, then C" and D' satisfy the twisted 
condition, i.e. we have 

Hb{c) X -ffb(D) = 0- 
Proof: The base matrix Hb{c) is a quasi-cyclic LDPC 
matrix which has the model matrix same as that of C with 
circulant matrices of size 1. By Proposition 13.11 C" and D' 
satisfy the twisted condition. 

Remember that two linear codes satisfy the twisted condi- 
tion if and only if these parity-check matrices Hi and H2 are 
orthogonal to each other i.e. 

Hix H^ = 0. 

Thus we have -ffb(c) x ^"[{0) — ^- " 

IV. NON-EXISTANCE RESULTS FOR A PAIR OF QC-LDPC 

Codes of Type IEEE802. 16e with Rate 1/2 and 2/3B 

The column (resp. row) weight distribution is defined as the 
distribution of Hamming weights of the columns (resp. the 
rows) for a given parity-check matrix. The weight distribution 
is one of the important parameters for the design of LDPC 
codes. In fact, the performance of LDPC codes with well- 
optimized weight distribution is very close to the asymptotic 
theoretical bounds [10]. Note that the weight distribution of the 
parity-check matrix of quasi-cyclic LDPC code with circulant 
matrices is the same as that of its base matrix. There are 
columns (resp. rows) of Hamming weight w in the base matrix 
if and only if there i^z columns (resp. rows) of Hamming 
weight w in associated parity-check matrix, where z is the 
size of the circulant matrices. 

In this section, we discuss our problem, which is a con- 
struction of a pair of irregular LDPC codes to derive a CSS 
code, under the following conditions: 

> Classical code C and D are quasi-cyclic LDPC codes 
with circulant matrices of size z. 

• The associated low-density parity-check matrices of C 
and D have the same "row-weight" distributions. 

• The weight of columns is more than or equal to 2. This 
condition arises from a standard decoding method, which 
is called a sum-product decoding. If one of a column- 
weight is less than 2, then the sum-product decoding does 
not work well [12]. 



In this paper, denote the condition above by (I). The condition 
(I) does not impose the Hmitation on the base matrices. 

In subsection i ^IV-AI we discuss the impossibility on con- 
struction of D under i?b(c) ■= -H^b(2/3B)- Similarly, in subsec- 
tion i llV-BI the impossibility on construction under H^c) '■= 
Hb(i/2) is discussed. 

A. Construction Impossibility on D under H^qj := i?{,(2/3B) 
Let C be a quasi-cyclic LDPC code with the base matrix 



b{C) 



6(2/35)- 



It is required that the base matrix of D satisfies: 

Hb{c) X H^{D) = 0- 
It implies that any row x G Fj'' of satisfies: 



Hb(c)x = 0. 



Define a matrix M as follows: 



M 



11111111 
01111111 
00111111 
00011111 
00001111 
00000111 
00000011 
11111110 



Then, for a vector x G F2'', we have 



MH, 



By an easy calculation, we obtain 



MHbt^c) = 



000000000000000010000000 
101010101010101001000000 
111111111111111100100000 
010101010101010100010000 
000000000000000000001000 
101010101010101000000100 
111111111111111100000010 
010101010101010100000001 



By considering the first and the fifth rows, we obtain a 
necessary condition 2:17 = X21 = for H^c)X = 0, where 

X = {xi,X2, ■ . ■,X24)- 

Thus the weight of 17th column and 21st column of -fffc(D) 
must be 0, in particular is less than 2. Therefore it is impossible 
to construct D which satisfies the condition (I) under ^^6(c) — 
Hb(2/3B)- By the result, we immediately have the 

Corollary 4.1: Let C be one of 19 quasi-cyclic LDPC codes 
which has the "base matrix" written as rate 1/2 codes in 
IEEE802.16e. Then there is no quasi-cyclic LDPC code D 
which satisfies (I). 

Proof: C is a particular example which satisfies our 
condition. ■ 



B. Construction Impossibility on D under Hb(C) '■— -^6(1/2) 
Let C be a quasi-cyclic LDPC code with its base matrix 



b{C) 



b(l/2) 



By the twisted condition, for the base matrix Hb(D) of a 
quasi-cyclic LDPC code D, the following equation is required: 



Hb{C) X -ffb(D) 



0. 



It implies that any row x E Fj"* of Hb[D) satisfies: 



b(C)^ 



0. 



The weight of each row of Hb(D) is 6 or 7. By the definition 
of a quasi-cyclic LDPC code of type IEEE802.16e, two non- 
zero entries of any row of Hb^o)^ except for the below, appear 
in adjacent position to each other. The weight of the right-side 



11 bits of a row of Hi 



b{D) 



is one or two and we denote the 



weight of the right-side by wt|i3_. 
Define three sets Xi, X2 and Xj,: 



Xi := {x e F^'*|wt(2;) = 6 or 7}, 



X2 := {x e 
^3 



,24 



Xi — Xi^i = 1 for some 14 < z < 23}, 



{x e Ff |wt|i3_(a;) = 1 or 2}. 



For simplicity, we denote X :— X 1(^X2 HX^. Then each row 
of Hb(D) belongs to the set X. 

We can determine the set X by using a personal com- 
puter. In fact, we verify X consists of six elements. The 
right-side 12-bits of each element of X have the form 
(0, 0, 0, 0, 0, 0, 0, 0, 1, 1, 0, 0). On the other hand, eleven kinds 
of the right-side 12-bits are required to satisfy (I). It shows that 
it is impossible to construct Hb[D) which satisfies (I) under 
the assumption Hb[c) ■= ^b(i/2)- 

By the result, we immediately have the 

Corollary 4.2: Let C be one of 19 quasi-cyclic LDPC codes 
which has the "base matrix" written as rate 2/3B codes in 
IEEE802.16e. Then there is no quasi-cyclic LDPC code D 
which satisfies (I). 

Proof: C is a particular example which satisfies our 
condition. ■ 

V. Conclusion 

In this paper, we discuss the conditions to construct a pair of 
quasi-cyclic LDPC codes, in particular of type IEEE802.16e, 
to derive a CSS code. The key of our discussion is to 
analyze the base matrices of quasi-cyclic LDPC codes. By our 
research, we find the impossibility to construct a pair of quasi- 
cyclic LDPC codes to derive a CSS code under conditions 
that (I) holds and one of the base matrices is the same as an 
IEEE802.16e LDPC base matrix with rate 1/2 or 2/3B. Note 
that our results are more general. The 32 quasi-cyclic LDPC 
codes, which have the model matrice written as rate 1/2 and 
rate 2/3B in IEEE802.16e, are typical examples of our results. 

We should remember that six base matrices, totally 96 quasi- 
cyclic LDPC codes, are discribed in the standardization form 
of IEEE802.16e. Hope has been left for the other four LDPC 
codes, totally 76 quasi-cyclic LDPC codes. We leave these 
posibilities as open problems. 
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